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ABSTRACT 


Standing sausage modes in flare loops are important for interpreting quasi- 
periodic pulsations (QPPs) in solar flare lightcurves. We propose an inversion 
scheme that consistently uses their periods P and damping times T to diagnose 
flare loop parameters. We derive a generic dispersion relation governing linear 
sausage waves in pressure-less straight tubes, for which the transverse density 
inhomogeneity takes place in a layer of arbitrary width / and is of arbitrary form. 
We find that P and 7 depend on the combination of [R/vai, L/R,1/R, pi/ pel, 
where R is the loop radius, L is the looplength, va; is the internal Alfvén speed, 
and ein, is the density contrast. For all the density profiles examined, P and 


EE 2015: 


T experience saturation when L/R > 1, yielding an inversion curve in the 
[R/vai,l/R, pi/ pe] space with a specific density profile when L/R is sufficiently 
large. When applied to a spatially unresolved QPP event, the scheme yields that 
R/vai is the best constrained, whereas L/R corresponds to the other extreme. 
For spatially resolved QPPs, while L/R >> 1 cannot be assumed beforehand, an 
inversion curve remains possible due to additional geometrical constraints. When 
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a spatially resolved QPP event involves another mode, as is the case for a recent 
event, the full set of [vai,/, pi/pe] can be inferred. We conclude that the proposed 
scheme provides a useful tool for magneto-seismologically exploiting QPPs. 


Subject headings: magnetohydrodynamics (MHD) — Sun: flares — Sun: corona 
— Sun: magnetic fields — waves 


1. INTRODUCTION 


The original ideas that laid the foundation for the field of solar magneto-seismology 
(SMS) were put forward in the 1970s (Uchida 1970, Rosenberg 1970, Zaitsev & Stepanov 1975, 
see also Roberts et al. 1984). However, this field flourished only after a rich variety of low- 
frequency Magnetohydrodynamic (MHD) waves and oscillations were identified with the 
advent of the TRACE, SOHO, Hinode, and SDO satellites (for recent reviews, see Banerjee 
et al. 2007, De Moortel & Nakariakov 2012, Mathioudakis et al. 2013; and also Ballester 
et al. 2007, Nakariakov & Erdélyi 2009, Erdélyi & Goossens 2011 for three recent topi- 
cal issues). It is also indispensable to refine the theoretical understanding of the collective 
wave modes supported by magnetized tubes, thereby enabling one to employ the measured 
wave properties to infer the solar atmospheric parameters that are difficult to measure di- 
rectly (e.g., Roberts 2000; Nakariakov & Verwichte 2005). Regarding its applications to the 
solar corona, SMS can offer such key information as the magnetic field strength in coronal 
loops (e.g., Nakariakov & Ofman 2001; Erdélyi & Taroyan 2008; Ofman & Wang 2008; White 
& Verwichte 2012) and above streamer stalks (Chen et al. 2010, 2011), the magnitude of field- 
aligned loop flows (Li et al. 2013; Chen et al. 2014), the temperature of loop plasmas (e.g., 
Marsh & Walsh 2009), the coronal effective adiabatic index (Van Doorsselaere et al. 2011), 
as well as the longitudinal (Verth & Erdélyi 2008; Andries et al. 2009; Luna-Cardozo et al. 
2012a) and transverse structuring (e.g., Arregui et al. 2007; Goossens et al. 2008; Yuan et al. 
2015). In addition, SMS applications with torsional Alfvén waves have proven invaluable 
in inferring the magnetic field structure at chromospheric heights (Jess et al. 2009; Fedun 
et al. 2011). Likewise, Luna-Cardozo et al. (2012a) demonstrated the potential of using 
longitudinal waves to infer the longitudinal variation of density and magnetic field strength 
in chromospheric waveguides. 


Magneto-seismological applications with standing kink modes (with azimuthal wavenum- 
ber m = 1) have been a common practice since their detection with TRACE (Aschwanden 
et al. 1999). Kink oscillations tend to experience substantial damping (e.g., Ofman & As- 
chwanden 2002; Ruderman & Erdélyi 2009; Verwichte et al. 2013, and references therein), 
which is usually interpreted in terms of resonant absorption (Ruderman & Roberts 2002, 
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Goossens et al. 2002, also Hollweg & Yang 1988, and the comprehensive review by Goossens 
et al. 2011). With this interpretation, Ruderman & Roberts (2002) and Goossens et al. 
(2002) suggested that the measured period P and damping time 7 can be used to infer 
the lengthscale l of the density inhomogeneity across coronal loops in units of loop radius 
R. For this purpose, the largely unknown transverse density distribution was shown to be 
important, since its formulation has a considerable impact on P and 7 (Soler et al. 2013, 
2014). 


While kink modes have attracted much attention, sausage modes (with m = 0) are 
equally important in SMS. In fact, sausage modes are even more important from the stand- 
point of solar atmospheric heating given their stronger compressibility and ubiquity in the 
lower solar atmosphere (Morton et al. 2011, 2012; Freij et al. 2014; Dorotovié et al. 2014; 
Grant et al. 2015; Moreels et al. 2015). In addition, sausage modes are important for inter- 
preting quasi-periodic pulsations (QPPs) in the lightcurves of solar flares (see Nakariakov 
& Melnikov 2009, for a recent review). Two distinct regimes are known to exist, depending 
on the axial wavenumber k along flare loops (Nakariakov & Verwichte 2005). The trapped 
regime results when k exceeds some critical value ke, where the energy of sausage modes 
is well confined to magnetic tubes. When k < ke, the leaky regime arises and sausage 
modes experience apparent temporal damping by radiating their energy into the surround- 
ing fluid (Spruit 1982; Cally 1986). It is known that ke depends sensitively on the density 
contrast between loops and their surroundings (e.g., Kopylova et al. 2007). In addition, both 
eigen-mode analyses (Kopylova et al. 2007; Vasheghani Farahani et al. 2014) and numerical 
simulations from an initial-value-problem perspective (Nakariakov et al. 2012; Chen et al. 
2015) indicated that the period P of sausage modes increases smoothly with decreasing k (or 
equivalently with increasing looplength L given that k = a/L for fundamental modes) until 
reaching some saturation value P, for sufficiently thin loops (R/L « 1). Likewise, identically 
infinite in the trapped regime for ideal MHD fluids, the attenuation time 7 decreases with 
decreasing k before experiencing saturation at Ts when R/L < 1. 


Magneto-seismological applications of sausage modes are possible due to their depen- 
dence on atmospheric parameters (Luna-Cardozo et al. 2012a,b). The practice based on the 
measured period and damping time can be illustrated by the study presented in Kopylova 
et al. (2007) where a step-function (top-hat) form was adopted for the transverse density 
distribution. The saturation values, P, and 7,/P,, for large density contrasts are approxi- 
mately 
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As an example, Kopylova et al. (2007) examined the QPP in the radio emissions reported 
in McLean & Sheridan (1973), where P and 7/P were found to be ~ 4.3 secs and ~ 10, 
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respectively. With the damping attributed to wave leakage, Eq. (1) then yields a density 
contrast p;/p. of ~ 100, and a transverse Alfvén transit time Rina of ~ 1.64 secs, provided 
that the flaring loop in question is sufficiently thin. However, the dependence of P and 
T on loop parameters is substantially more involved if one goes a step closer to reality by 
replacing a step-function density profile with a smooth one. Even for thin loops where neither 
P nor T depends on looplength, both P, (Nakariakov et al. 2012) and 7, (Chen et al. 2015) 
may be sensitive to the steepness, or equivalently the lengthscale, of the transverse density 
distribution. In mathematical terms, this means that for a given density profile, 
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When only P, and 7, are known as is the case for measurements without imaging capabilities, 
the appearance of L/R no longer allows a unique pair of [p;/p., R/vai] to be deduced. Despite 
this, one can still constrain the combination [R/vai, pi/Pe,!/R] by developing a scheme in 
much the same way that kink modes were employed (Arregui et al. 2007; Goossens et al. 
2008; Soler et al. 2014), the only difference being that the transverse Alfvén time Rina 
replaces the longitudinal one Lia, 


The present study aims to develop the aforementioned scheme employing measured 
periods and damping times of standing sausage modes. An essential ingredient will be 
to establish the functions F and G in Eq. (2). To this end, we will derive an analytical 
dispersion relation (DR) governing linear sausage waves hosted by magnetized tubes with 
a rather general transverse density distribution. The only requirement here is that this 
density distribution can be decomposed into a uniform cord, a uniform external medium, 
and a transition layer connecting the two. However, the density distribution in the transition 
layer is allowed to be arbitrary, thereby making the DR applicable to a rich variety of density 
profiles. We note that this kind of density profiles has been extensively adopted in kink mode 
studies (e.g., Soler et al. 2014, and references therein). We further note that developing an 
analytical DR is important in its own right. Apart from the step-function profile (e.g., Spruit 
1982; Cally 1986), analytical DRs in the cylindrical case are available only for a limited set 
of density profiles (Edwin & Roberts 1986, 1988; Lopin & Nagorny 2014). 


This manuscript is organized as follows. Section 2 presents the derivation of the DR and 
our solution method. A parameter study is presented in Sect. 3 to examine how the period 
and damping time of sausage modes depend on tube parameters, thereby establishing our 
numerical scheme for inverting measurements of spatially unresolved QPPs. An extension 
to spatially resolved QPPs is then given in Sect. 4. Finally, Sect. 5 closes this manuscript 
with our summary and some concluding remarks. 
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2. MATHEMATICAL FORMULATION 
2.1. Description for the Equilibrium Tube 


We consider sausage waves in a structured corona modeled by a density-enhanced cylin- 
der with radius œR aligned with a uniform magnetic field B = BZ, where a cylindrical coor- 
dinate system (r,0,z) is adopted. The equilibrium density is assumed to be a function of r 
only and of the form 


Pi, 0<r<r = R-1/2, 
pr) = Gell Tri <r <růe= R+1/2, (3) 
Pes T> Tg: 


The profile between Ir, re] is such that the equilibrium density p decreases continuously from 
the internal value p; to the external one pe. The thickness of this transition layer, denoted 
by J, is bounded by 0 and 2R. The former represents the steepest profile of a step-function 
form, whereas the latter corresponds to the least steep case. 


While our analysis is valid for arbitrary prescriptions of pẹ, a number of choices have to 
be made to evaluate quantitatively the effects of equilibrium density profiles. To this end, 
we select the following profiles, 
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The profiles labeled linear, parabolic and sine have been examined in substantial detail in 
the context of standing kink modes (Soler et al. 2013, 2014). An additional profile, labeled 
inverse-parabolic, is added to make the list more comprehensive in that it naturally com- 
plements the parabolic one. Figure 1 illustrates the r-dependence of the chosen equilibrium 
density profiles. For illustration purposes, pi/pe is chosen to be 100, and //R is chosen to be 
unity. 
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2.2. Solutions for Radial Lagrangian Displacement and Total Pressure 
Perturbation 


Appropriate for the solar corona, we work in the framework of cold (zero-3) MHD, in 
which case sausage waves do not perturb the z-component of the plasma velocity. Let dv, 
denote the radial velocity perturbation, and let ôb, and 6b, denote the radial and longitudinal 
components of the perturbed magnetic field db, respectively. The perturbed total pressure 
is then Ana = B - 6b/4a = Bdb,/4a given the absence of thermal pressure in the zero-( 
limit. Fourier-decomposing any perturbed value df (r, z; t) as 


pts, z;t) = Re { F(r) exp [i (wt — ul l (5) 


one finds from linearized, ideal, cold MHD equations that 
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where the prime’ = d/dr. In addition, € = i#,/w is the Fourier amplitude of the radial La- 
grangian displacement, and va(r) = B/,/47p(r) is the Alfvén speed. The Fourier amplitude 
of the perturbed total pressure is 
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With azimuthal wavenumber m being 0, the equations governing linear sausage waves 
are free of singularities, making our derivation simpler than in kink mode studies where a 
treatment of singularity is necessary (Soler et al. 2013). To be specific, the solutions to 
Eq. (6) in the transition layer can be expressed as a regular series expansion in x = r — R. 
Let bmn and Exo denote two independent solutions, 
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Without loss of generality, one may choose Joo, ai] = [R, 0] and [bo, bi] = [0, 1]. Expanding 
the equilibrium density p about x = 0 as well, one finds that 
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Plugging Eq. (8) into Eq. (6) with the change of independent variable from r to x, and then 
ordering the terms according to powers in x, one finds that 
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where 7 = 47/B? and x represents either a or b. 


With Eq. (8) at hand, the solution to Eq. (6) can be expressed as 


Ai Jy (ur), O<r<ri, 
Entr) = Auf a(x) ge Aaf 2) ry < r < Te; (12) 
AH (uer), F Z Te 


where Ai, Ae, A; and Ag are arbitrary constants, and J, and HY are the n-th-order Bessel 
and Hankel functions of the first kind, respectively (here n = 1). In addition, ui, = 
w*/Viie KI with vie = B?/(4mpie). As discussed in Cally (1986), requiring that —1/2 < 
arg Hi, arg He < 7/2 does not exclude any additional independent solution. Furthermore, ex- 


pressing the external solution in terms of H UI 


permits a unified examination of both trapped 
and leaky waves. Indeed, the trapped regime arises when arg He = 7/2, in which case one 
finds that H® (uer) = —(2/m)Ky(ar) with a = e/i being real and positive (see discussions 
on page 281 in Cally 1986). Now the Fourier amplitude for the total pressure perturbation 


can be evaluated with Eq. (7), the results being 
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The expression for Prot for the ranges r < ri and r > rẹ can also be independently verified 
by using an alternative relation between €, and Prot, 
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2.3. Dispersion Relation of Sausage Waves 


The dispersion relation (DR) governing linear sausage waves can be derived by requiring 
that both €, and pi. be continuous at the interfaces r = ri and r = re. This leads to 
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where x; = —l/2 and x, = 1/2. Eliminating A; (Ae) by dividing the third (fourth) by the 
first (second) equation, one finds that 


AA + AgAe = 0, 
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Evidently, for Eq. (14) to allow non-trivial solutions of (Av, A2], one needs to require that 
Ay Ag —_ AgAg3 = 0, (17) 


which is the DR governing sausage waves in nonuniform loops. 
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Does Eq. (17), the DR valid for arbitrary //R, recover the well-known result for the 
step-function profile when //R — 0? Retaining only terms to the 0-th order in //R and 
noting that ri ~ re © R, one finds that A, (n = 1,--- ,4) simplify to 


Ay = (1 — Xj)ao + Ray, 
A = (1 — Xi)bo + Rb, 
A = (1 — Elan + Ray, 
Ag = (1 — Xe)bo + Rb. 


Substituting these expressions into Eq. (17), one finds that 
(Xi = Xe)(aıbo = oof) = 0. 


This leads to X; = Xe, given that abọ — gan is not allowed to be zero for ra and D to 
be independent. In other words (see Eq. (16)) 
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which is the DR for equilibrium density profiles of a step-function form (e.g., Cally 1986). 


While the DR is equally applicable to propagating waves, we will focus on standing 
modes for which the axial wavenumber k is real, while the angular frequency w is allowed to 
be complex-valued. Furthermore, let us focus on fundamental standing modes supported by 
magnetized tubes of length L. In this case, another measure to validate the DR, independent 
of the eigen-value problem approach, is to employ the linearized, time-dependent, cold MHD 
equations to derive an equation governing the transverse velocity perturbation v,(r, z,t). In 
view of the line-tying boundary conditions at the loop ends z = 0 and z = L, one may 
express v,(r,z,t) as v(r,t)sin(kz) with k = zi, yielding (Nakariakov et al. 2012; Chen 


et al. 2015) 
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When supplemented with appropriate boundary and initial conditions, the signal of v,(r, t) 
at some arbitrarily chosen distance from the tube axis can be followed. As demonstrated 
in Nakariakov et al. (2012); Chen et al. (2015), after a transitory phase this signal evolves 
into a harmonic (decaying harmonic) form when k is larger (less) than some critical value, 
corresponding to the well-known trapped (leaky) regime. Numerically fitting the signal 
with a sinusoidal (exponentially decaying sinusoidal) function then yields the period P (P 
together with the damping time 7) for trapped (leaky) modes. We also adopt this approach 
and compare the derived values for P and r with what is found by solving the DR for 
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complex-valued w at given real k. As will be shown in Fig. 2, the two sets of independently 
derived values for [P,7] agree remarkably well. At this point, it should be remarked that 
once a choice for pẹ is made, P and T depend only on the combination of parameters 
[pi/ Pe, !/R, L/R] when they are measured in units of the internal Alfvén transit time R/vai. 
Here we have used L/R = 7/(kR) in place of the dimensionless axial wavenumber kR. 


In general, the DR (Eq. (17)) is not analytically tractable and is solved numerically 
for a given pẹ profile and some given combination of [pi/pe,l/R, L/R]. To do so requires 
the infinite series expansion in Eq.(8) to be truncated by retaining the terms with n up toa 
certain N. A value of N = 101 is chosen for all the numerical results to be presented, and we 
have made sure that using an even larger N does not introduce any appreciable difference. 
In addition, we focus only on the lowest order modes, namely those with the simplest radial 
structure in the eigen-functions. When verifying these results with the computations from 
an initial-value-problem perspective, we choose not to use a too localized initial perturbation 
u(r,t = 0) (see Eq.(19)), otherwise higher order modes are introduced to contaminate the 
v(r,t) signals. 


3. NUMERICAL RESULTS AND THE INVERSION SCHEME 


Let us start with an examination of how the gu profile impacts the dispersive properties 
of standing sausage modes. Figure 2 shows the dependence on the length-to-radius ratio 
L/R of the period P and damping time 7 for different choices of pẹ(r) as labeled. For 
illustration purposes, we choose the density contrast p;/pe to be 100, and choose L/R to 
be unity. The black line in Fig. 2a, which represents 2L/VAe, separates the trapped (to 
its left, where 7 is identically infinite) from leaky (to its right) regimes. The curves are 
found by solving the analytical DR (Eq. (17)), whereas the open circles are found by solving 
the corresponding time-dependent equation (see Eq. (19) and the associated description). 
Evidently, the periods P and damping times 7 obtained from the two independent approaches 
agree with each other remarkably well. Figure 2 indicates that the overall tendency for P (rT) 
to increase (decrease) with L/R is seen for all the equilibrium density profiles considered. 
In particular, regardless of the profiles, both P and 7 tend to some asymptotic values at 
large L/R. However, the choice of equilibrium density profiles has a considerable influence 
on the specific values for P and 7. This is particularly true if one compares the results for 
the parabolic and inverse-parabolic profiles, given by the green and blue curves, respectively. 
Furthermore, while the periods P for the linear and sine profiles differ little for the chosen 
L/R, the damping times 7 show a stronger profile dependence. This signifies the importance 
of using P and 7 in a consistent manner when one attempts to deduce how the transverse 
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equilibrium density is structured. 


The effects of equilibrium density profile can be better brought out by capitalizing on 
the fact that both P and T experience saturation for sufficiently thin loops. Let P, and 
Ts denote the saturation values. Figure 3 presents how P, and 7,/P; depend on L/R, the 
width of the transition layer in units of loop radius. The density contrast is chosen to be 
100, and different choices of gu are represented by the curves in different colors as labeled 
in Fig. 3b. In addition, in place of Ts, the ratio 7,/P, is plotted since it is a better measure 
of signal quality. One sees from Figs. 3a and 3b that the curves converge at L/R — 0 as 
expected given that the DR (Eq. (17)) simplifies to Eq. (18) pertinent to a top-hat density 
distribution. Figure 3a indicates that the //R-dependence of P, critically depends on how 
Ptr is described. For the parabolic profile, P, increases monotonically with //R, whereas the 
opposite trend is found for the inverse-parabolic profile. When it comes to the linear and 
sine profiles, Fig. 3a shows that the //R-dependence of P, is not as strong, and the difference 
between the two profiles is discernible only when //R = 1. Moving on to Fig. 3b, one sees 
that the ratio 7,/P; decreases monotonically with //R for all the profiles, meaning that wave 
leakage plays an increasingly important role in attenuating sausage modes when the loop 
becomes more diffuse. Reinforcing the impression from Fig. 2b, one sees that relative to P,, 
T;,/P; better discriminates the equilibrium density profiles. 


So far we have fixed the density contrast pi/pe at 100. One naturally asks what happens 
if p;/pe is varied? Figure 4 presents the distribution of P, (the left column) and 7,/P, (right) 
in the [pi/pe,l/R] plane. Each row represents one of the four density profiles as labeled. 
Besides, the red curve represents where 7,/P; = 10, the value for the QPP event reported 
in McLean & Sheridan (1973). Examine the left column first. One sees from Figs. 4b and 
4c that in the parameter range examined, P, tends to increase (decrease) with //R at any 
given pi/pe when the parabolic (inverse-parabolic) profile is chosen. In contrast, Figs. 4a and 
4d indicate that P, for the linear and sine profiles shows a nonmonotonical dependence on 
L/R, even though this variation is hardly discernible. Now consider the right column, from 
which one can see that regardless of the profiles, 7,/P, decreases with increasing L/R for all 
the density contrasts examined. In addition, the dependence of 7,/P, on L/R is the strongest 
for the inverse-parabolic profile, and the least strong for the parabolic one. The dependence 
for the linear and sine profiles lies in between, with the dependence in the linear case being 
slightly stronger. 


Conceptually, Fig. 4 can be used to invert the measured values of the period and damping 
time of sausage modes, provided that the loops hosting these oscillations are sufficiently thin. 
Consider the QPP event reported in McLean & Sheridan (1973) as an example, for which 
P, = 4.3 secs and 7,/P, = 10. It then follows that for a given density profile, any point 
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along the corresponding red curve in Fig. 4 can equally reproduce the measured 7,/P,. After 
reading any pair of [p;/p.,1/R], one can read from the left column the corresponding value 
for P, in units of Riva, With P, known, one can then deduce Rina, In practice, however, 
constructing such a contour plot is not necessary, and one may simply consider the following 
3-step inversion scheme. First, one starts with the dispersion relation for a step-function 
density profile (L/R = 0), Eq. (18), to find the value for the density contrast p;/p. such that 
T| P; agrees with the measured value. Second, with the p;/p. value for a smaller //R as a 
good guess, one can then solve Eq. (17) to find a new pi/pe that yields the measured 7,/P, 
by increasing //R from 0 to 2 consecutively. Third, with pi/pe found for all possible //R, one 
can solve Eq. (17) for a given pair of [pi/pe,l/R], yielding a value for P,/(R/va;). Finding 
the transverse Alfvén transit time R/va; is then straightforward since P, is known. 


The product of the inversion scheme is an inversion curve in the three-dimensional (3D) 
space formed by Rina, L/R and p;/pe. Figure 5 presents such curves (the solid lines) and 
their projections onto various planes (dashed) for the examined density profiles, pertinent 
to the event reported in McLean & Sheridan (1973). To help digest this figure, a number 
of points are read from the curves and presented in Table 1. One sees that among the 
parameters forming this 3D space, R/va; is the best constrained. The biggest (smallest) 
value, 2.13 (1.18) secs, is found for the inverse-parabolic (parabolic) profile when //R > 2. 
In other words, the biggest value exceeds the smallest one by only 79.7%. As to the density 
contrast, the biggest value (251.4 found for the inverse-parabolic profile when (UR — 2) is 
larger than the smallest one (88.1 when //R — 0) by 185%. The least constrained parameter 
is L/R, with any value in the allowed range from 0 to 2 being possible. 


4. FURTHER DEVELOPMENT OF THE INVERSION SCHEME 


Before proceeding, let us first recap the key points in the scheme for inverting the 
measured period P and damping time T of sausage modes. From the outset, we have assumed 
that only P and 7 are known, whereas we have no information on either the geometric 
parameters (R, L, and l) or the physical parameters (vai, pi, and pe). On top of that, 
the specific density profile is also assumed to be unknown. This happens when one has only 
spatially unresolved observations (see e.g., the majority of the events compiled in Aschwanden 
et al. 2004, Table 1). In this case, the dispersion relation, Eq. (17), suggests that in general 
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the periods and damping times can be formally expressed as 
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With only two measured values available, any point on a two-dimensional (2D) surface 
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in the four-dimensional (4D) space formed by [R/vai, L/R,l/R, ein, can reproduce the 
measurements, even if one is allowed to prescribe a density profile. If the sausage modes are 
in the trapped regime in the sense that the QPPs do not show temporal damping, then the 
situation becomes even less desired since now the restriction from 7 is no longer available. 
This complexity can be alleviated if the condition L/R >> 1 holds for the flare loops in 
question since L/R no longer appears such that Eq. (2) is restored. One then finds a curve 
in the [R/vai,1/R, pi/pe] space as shown by Fig. 5. 


The situation improves if the QPP events are spatially resolved, since the looplength 
(L) and the outer interface of the loops (re = R + 1/2) can be considered known. In this 
case L/R and L/R are no longer independent but are related by 


L L = l (21) 
R r 2bRI" 
where we have used the relation 
Kes 
R= ——. 22 
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Now the situation is similar to what Eq. (2) implies: if the sausage mode is a trapped (leaky) 
one, then a 2D surface (1D curve) can be deduced in the [R/vai, L/R, p;/pe] space. It is just 
that now for a point on this 2D surface or along the 1D curve, one can further deduce va; in 
view of Eq. (22). 


Something fascinating happens if one observes a spatially resolved QPP event hosting 
more than just a sausage mode. For illustration purposes, we consider the situation where a 
temporally decaying kink mode is involved, and its damping can be attributed to resonant 
absorption. In this case, the period and damping time for the kink mode can be formally 
expressed by 


= {4 


Note that while L/R in principle can be incorporated into Eq. (23), in reality there is no 
need to do so since the corrections to Pkink and Tkink due to finite L/R are of the order 
(R/L)? (Van Doorsselaere et al. 2004; Goossens et al. 2008). Even for relatively thick flare 
loops, R/L is of the order 0.1 and these corrections amount to only a few percent. Note 
further that Fkink and Gkink have been extensively studied, and a graphical representation 
can be found in Fig. 1 of Soler et al. (2014). With Fink, Grink, Faus, and Gaus known, it 
is then possible to fully constrain the unknowns |I, R, vai, pi/Pel, if one assumes a density 
profile. In fact, the measured values of |Pkink, Tkink, Paus: Tsaus] are more than sufficient: 
l and R are not independent but are related through Eq. (22). This suggests that only 
three expressions contained in Eqs. (20) and (23) are needed. In practice, we consider the 
expression for the kink mode period as the redundant one. Now the inversion procedure is 
rather straightforward. In view of the relation (21), the two expressions for the damping- 
time-to-period ratio in Eqs. (20) and Eqs. (23) contain only two unknowns, namely //R 
and p;/pe. With both Tyaus/Peaus ANd Tkink/Pkink available, one can deduce a unique pair 
of [1/R,p;/p-]|. The loop radius R then follows from Eq. (22), which then enables one to 
evaluate va; with the first expression in Eq. (20). Finally, as a safety check, one can proceed 
to evaluate, with the first expression in Eq. (23), the theoretically expected kink mode period 
Prink,theory: The deviation of Pyink,theory from the measured one then allows to say a few words 
on how safe it is to identify the oscillation signals with some particular modes. 


As an illustration of the aforementioned inversion procedure, let us consider the QPP 
event in microwave emissions measured with the Nobeyama Radioheliograph (NoRH) on 
14 May 2013 (Kolotkov et al. 2015). Lucky enough, it is likely that this event contains 
a fundamental kink mode with Pkink = 100 secs and Tkink/Pkink = 2-5, in addition to a 
fundamental sausage mode with Paus = 15 secs and Tgaus/Psaus = 6. Assuming that the 
apparent width measured therein corresponds to 2re, one finds that re = 4 x 10° km and 
L = 4 x 104 km, meaning that L/re = 10. Let us assume that wave leakage is responsible 
for damping the sausage mode, and resonant absorption is responsible for damping the kink 
one. Furthermore, let us assume that the sine profile best describes the equilibrium density 
distribution. In this case, the analytical expressions obtained in the thin-tube-thin-boundary 
approximation are accurate to within ~ 25% (e.g., Van Doorsselaere et al. 2004; Soler et al. 
2014). This enables us to illustrate our inversion procedure without resorting to a fully 
numerical solver to establish Fkink and Gink. Now the formal expressions given by Eq. (23) 
can be replaced with (e.g., Goossens et al. 2008) 


kink =~ = SSES 
VAi pil Pe (24) 
Tkink _ 2 pi/Pe+1 1 


est 8 T pif Pe — 11/R’ 
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Following the outlined inversion procedure, we find that L/R = 0.272 and p;/pe = 29.8 as 
constrained by the ratios T/P. It then follows that R = 3.52 x 10° km, from which one can 
deduce that va; = 623 km s `. Finally, Pyink theory is found to be 92.2 secs, which agrees 
with the measured value of 100 secs to within 8%. This safety check lends support to the 
interpretation of the two modes in terms of fundamental kink and sausage modes as done 
by Kolotkov et al. (2015). Besides, the deduced Alfvén speed va; and density contrast pi/pe 
both seem reasonable. On top of that, there is no need for one to worry too much about 
the accuracy of Eq. (24) describing the kink mode: at the deduced [I/R, p;/pe], this equation 
yields P and 7/P to an accuracy better than 1% and ~ 6%, respectively (see Fig. 1 in Soler 
et al. 2014). 


5. CONCLUSIONS 


How plasma density is structured across various magnetic structures in the solar corona 
remains largely unknown. It has been a common practice to deduce this key information by 
employing magneto-seismological techniques that invert the measured period P and damping 
time 7 of standing kink modes collectively supported by a magnetic structure (e.g., Goossens 
et al. 2008; Soler et al. 2014). In contrast, while quasi-periodic pulsations (QPPs) in the 
lightcurves of solar flares are often attributed to standing sausage modes in flare loops and 
therefore can also offer the associated period and damping time, a scheme is missing for 
inverting these two measurements to deduce the information on the density distribution 
transverse to flare loops. The primary aim of this study has been to construct such a 
scheme. To this end, we worked in the framework of cold (zero-3) MHD and modeled flare 
loops as straight cylinders with a transverse density profile characterized by a transition 
layer sandwiched between a uniform cord and a uniform external medium. An analytical 
dispersion relation (DR) governing linear sausage waves, Eq. (17), was derived by solving the 
perturbation equations in terms of a regular series expansion in the transition layer. This 
DR, valid for arbitrary choices of the density profile in the transition layer, then enabled us 
to examine the effects of density structuring on the periods and damping times of sausage 
modes, thereby facilitating the construction of the inversion scheme. 


In general, we found that P and T of sausage modes depend on a combination of param- 
eters [R/vai, L/R,1/R, p;/pe| as formally expressed by Eq. (20), where the functions Fyaus 
and Gaus are a product of the DR. Here L and R denote the looplength and loop radius, 
respectively. Furthermore, / is the width of the transition layer, va; is the Alfvén speed in 
the cord, and pi/pe is the density contrast between the loop and its surrounding fluid. We 
showed that for the density profiles examined, both P and 7 experience saturation for suffi- 
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ciently large L/R when the rest of the four parameters are fixed. The choice of the transverse 
density profile was found to have a considerable influence on P and 7, their dependence on 
L/R in particular. 


Our inversion scheme can find applications to both spatially unresolved and resolved 
QPP events. For spatially unresolved ones, we showed that the best one can do is to deduce 
a 1D curve in the [R/vai, L/R, pi/pe] space for a prescribed density profile. This happens if 
the QPPs in question experience temporal damping, and if the flare loops hosting them can 
be assumed to be sufficiently thin. When applied to a QPP event reported by McLean & 
Sheridan (1973), this inversion technique indicated that the transverse Alfvén transit time 
R/vaj is the best constrained, varying by a factor of 80% if the uncertainties in specifying the 
density profile are taken into account. The density contrast pi/pe is less well constrained, 
with the largest deduced value exceeding the smallest one by a factor of 1.85. The least 
constrained is the transverse density lengthscale in units of loop radius (J/R), any value in 
the allowed range (0,2) can be equally possible to reproduce the measurements. 


When it comes to spatially resolved events, the geometric parameters L and R+ 1/2 are 
additional constraints, on top of the measured values for P and T. Even though in this case 
one cannot assume L/R > 1 a priori, it is possible to deduce a 1D curve in the Ia, l, pi/pe] 
space for a chosen density profile since R is expressible in terms of L/R. If a spatially resolved 
QPP event comprises more than just one sausage mode, then it is possible to deduce the 
full information on |l, R, vai, pi/pe]. In this case, the inversion problem becomes an over- 
determined one. We illustrated this fascinating application with the QPP event reported 
by Kolotkov et al. (2015) where a fundamental kink and a fundamental sausage mode were 
suggested to co-exist and both experience temporal damping. Attributing the temporal 
damping of the kink mode to resonant absorption, and that of the sausage mode to wave 
leakage, we were able to deduce the full set of [l, R, vai, pi/pe] by using Eqs. (20) and (24) 
and assuming a sine profile for the density distribution. One redundant equation, here taken 
as the expression for the kink mode period, can then allow a safety check on, say, whether 
it is reasonable to interpret the signals in the QPP event as the aforementioned modes. For 
this particular event, our results demonstrated that the interpretation provided in Kolotkov 
et al. (2015) is reasonable, not only because the deduced parameters seem realistic, but also 
because the theoretical prediction for the kink mode period agrees with the measured one 
remarkably well. 


Our scheme nonetheless has a number of limitations. First, we have assumed that 
the temporal damping of sausage modes is due to wave leakage, an ideal MHD process. 
In reality, non-ideal mechanisms like electron heat conduction and ion viscosity can provide 
additional channels for damping sausage modes. While these non-ideal processes were shown 
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by Kopylova et al. (2007) to be unlikely the cause for the temporal damping in the QPP 
event reported by McLean & Sheridan (1973), their importance has yet to be assessed for 
the event reported by Kolotkov et al. (2015). Second, working in the cold MHD, we have not 
taken into account the possible effects due to finite plasma beta, which may be of the order 
unity in flare loops. However, the corrections due to finite beta seem marginal (Inglis et al. 
2009). Third, the longitudinal variation in neither the plasma density nor the magnetic field 
strength has been considered, even though the corrections due to this variation are unlikely 
to be significant (Pascoe et al. 2009). Fourth, the density inhomogeneity in flare loops was 
assumed to be in a monolithic form, whereas in reality these loops may be multi-stranded. 
While the fine structuring in the form of randomly distributed concentric shells is found to 
have a far less significant influence than the monolithic component of the density distribution 
(Chen et al. 2015, see also Pascoe et al. 2007), there is a need to rigorously assess the effects 
due to fine structuring in the form of randomly distributed strands. 


We thank the referee for the constructive comments, which helped improve this manuscript 
substantially. This work is supported by the National Natural Science Foundation of China 
(41174154, 41274176, and 41474149), and by the Provincial Natural Science Foundation of 
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Fig. 1.— Transverse equilibrium density profiles as a function of r. The profiles differ only 


in how they are described in a transition layer sandwiched between the internal (with a 
uniform density e) and external (with a uniform density pe) parts. The transition layer is 
of width I, and is located between r; = R—1/2 and re = R+1/2, with R being the cylinder 
radius. Four different choices of the density profiles in the transition layer are adopted as 
labeled, and are given by Eq. (4). For illustration purposes, l is chosen to be R, and p;/p¢ is 
chosen to be 100. 
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Fig. 2.— (a) Period P and (b) damping time 7 as functions of loop length L. Four 
different choices of density profiles are examined as labeled. The black line in (a) represents 
P = 2L/vae and separates the trapped (to its left) from leaky (right) regimes. The open 
circles represent the values for P and 7 obtained by solving Eq. (19) from an initial-value- 
problem perspective, which is independent from the eigen-value-problem approach presented 
in the text. Here the width of the transition layer | = R, and the density contrast p;/p. = 100. 
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Fig. 3.— Saturation values for (a) period P, and (b) damping time 7, as functions of 


the width of the transition layer l. These saturation values are attained for sufficiently thin 
loops. Four different choices of the density profiles are examined as labeled. Here the density 
contrast pi/Pe = 100. 
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Fig. 4— Contour plots in the [//R, pi/pe]| space for the saturation values for period P, (the 
left column) and damping-time-to-period ratio Ts/P, (right). These saturation values are 
attained for sufficiently thin loops. Each row represents one of the four different choices of 
the density profiles as labeled. The red curve in each panel represents where 7,/P; = 10, 
corresponding to the event reported in McLean & Sheridan (1973) 
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Fig. 5.— Inversion curves (the solid lines) together with their projections (dashed) in the 
three-dimensional space formed by [R/vai,!/R, pi/pe|. Four choices for the density profiles 
are examined and given in different panels as labeled. All points along an inversion curve are 
equally compatible with the quasi-periodic-pulsation event reported in McLean & Sheridan 
(1973) where the oscillation period is 4.3 secs, and the damping-time-to-period ratio is 10. 
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linear parabolic inverse-parabolic sine 


ae min, R/vmĘ (sec) min, R/vai (sec) | pin, R/vai (sec) | pi/pe R/vai (sec) 
0.01 88.1 1.62 
0.2 89.1 1.62 
0.4 91.9 1.62 
0.6 96.4 1.63 
0.8 102.4 1.63 
1.0 109.8 1.62 
Ge 118.2 1.62 
1.4 127.4 1.61 
1.6 136.9 1.6 
1.8 146.4 1.58 
1.99 155.4 1.56 


Table 1: The inverted values for the transverse density length scale in units of loop radius 
(L/R), density contrast p;/p. and the transverse Alfvén transit time R/va;. This inversion is 
made for the quasi-periodic-pulsation event reported in McLean & Sheridan (1973), assuming 
that the associated flare loop is sufficiently thin. For this event, the oscillation period is 4.3 
secs, and the damping-time-to-period ratio is 10. 


